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ABSTRACT 

Precipitation generates small-scale turbulent air flows the energy of which ultimately dissi- 
pates to heat. The power of this process has previously been estimated to be around 2-4 
W m -2 in the tropics: a value comparable in magnitude to the dynamic power of the global 
circulation. Here we suggest that this previous power estimate is approximately double the 
true figure. Our result reflects a revised evaluation of the mean precipitation path length 
Hp. We investigate the dependence of Hp on surface temperature, relative humidity, tem- 
perature lapse rate and degree of condensation in the ascending air. We find that the degree 
of condensation, defined as the relative change of the saturated water vapor mixing ratio in 
the region of condensation, is a major factor determining Hp. We estimate from theory that 
the mean large-scale rate of frictional dissipation associated with total precipitation in the 
tropics lies between 1 and 2 W m -2 and show that our estimate is supported by empirical 
evidence. We show that under terrestrial conditions frictional dissipation constitutes a mi- 
nor fraction of the dynamic power of condensation-induced atmospheric circulation, which 
is estimated to be at least 2.5 times larger. However, because H P increases with surface 
temperature T s , the rate of frictional dissipation would exceed that of condensation- induced 
dynamics, and thus block major circulation, at T s > 320 K in a moist adiabatic atmosphere. 
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1. Introduction 



Understanding the physics of a moist at mosphere and capturing it in theoretical concepts 



is a major challenge for climate science (ISchneider 



2006 



Schiermeier 



2010 



Cotton et al. 



20111 ). Among the complications introduced by water vapour are the various influences of 



precipitation on atmospheric motion. One specific aspect is that precipitation generates 
small-scale air turbulence around the falling condensate particles. The energy for this tur- 
bulent air motion derives from the potential energy of the rain drops or ice particles (i.e. 
" hydrometeors" ) in the gravitational field of Earth and is ultimately dissipated to heat. If 
this energy were not converted to turbulent kinetic energy of the air, the hydrometeors would 
continue to accelerate as they fall. But air exerts a drag force that prevents this acceleration. 
This force grows with increasing size of the hydrometeor and its velocity W relative to the 
surrounding air. Thus, as the condensate particle is accelerated by gravity, this opposing 
force grows until it equals the weight of the particle. Acceleration ceases at this resulting 
terminal velocity Wt- Since hydrometeors fall at near their terminal velocities for most of the 
duration of their falls, the mean drag force acting on them over this period is approximately 
equal to their weight. 

Consider a column of moist mixture of dry air and water vapor which we de- 

note here using subscripts d and v respectively, standing for dry air and water vapor. In 



hydrostatic equilibrium, 



dp 

oz 



where p = pd + p v is air pressure, p = pd + p v is air density, g is the acceleration of 
gravity. There is no available potential energy in such a column. We now cool this column 



in such a manner that some water vapor condenses. Two types of potential energy have now 
become available: the potential energy of droplets in the gravitational field and the potential 
energy of any non-equilibrium air pressure gradient that may have formed upon condensation. 
We recently proposed that the release of the second type of potential energy, i.e. that 
associated with the non-equilibrium gradient of saturated water vapor, is a major driver 



of atmospheric 


circul 


Gorshkov et al. 


2012) 



ation on Earth (IMakarieva and Gorshkov! 120071 ; iMakarieva et al.l 12010 



20121 ). Since precipitation and, hence, frictional dissipation of the potential 
energy of hydrometeors, always accompany condensation, it is important to estimate and 
contrast the power of the two processes in order to understand their relative influences. 
In this paper we first examine how the power D of precipitation-related frictional dissipa- 



tion can be estimated from 
those of 



Pauluis et al 



aasic atmospheric parameters. We then compare our results to 



( 120001 ) and explain why our estimates are m ore consistent with both 



Pauluis and Diasl ( 12012( 1 



theory and data. We also discuss why the recent estimate of D by 
from satellite-derived tropical rain rates does not constrain the true value more accurately 
than our theoretical analysis. We show that D grows with increasing surface temperature 
and estimate the critical temperature when the power of frictional dissipation equals the 
power of condensation-induced dynamics, such that the latter ceases. 
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2. Basic formulae 

The power of dissipation of energy D (W m -2 ) per unit area associated with small-scale 
turbulence around hydrometeors can be written as 

poo poo 

D= / WF c dz= / Wtgpcdz. (2) 
Jo Jo 

Here W = w c — w > is the mean velocity of hydrometeors relative to the air, w c and w 
are the vertical velocities of condensate particles and air relative to the Earth's surface, F c 
is the turbulent drag force per unit air volume exerted by air on hydrometeors, p c = N c m is 
condensate density, N c is the number of hydrometeors per unit volume and m is their mean 
mass. The second equality in (j2J) takes into account our assumption that hydrometeors are 
falling at their terminal velocity Wt, such that the drag force acting on a droplet is equal to 
its weight. For a constant Wt Eq. (J2]) simply represents the product of Wtg and the total 
amount of condensate in the atmospheric column. 

Equation (J2J) is not suited for a theoretical analysis. The distribution of W t , which 
depends strongly on particle size, is poorly known. The amount of condensate in the at- 
mosphere varies greatly in time and space, from a few kilograms per square meter in se- 
vere storms to le s s than a hundred grams per square meter under normal conditions (e.g., 



Jiang et al. 



2008 



Wood et al. 



20021 ). Nonetheless, as we shall now show, it is possible to 
modify Eq. (j2J) to exclude these uncertain parameters. 

Let us first consider the case when the terminal velocity of hydrometeors is much larger 
than air velocity, Wt 3> \w\. Neglecting for now re-evaporation of condensate in the column 
we assume that all condensate that has formed in the course of the ascent of moist air 
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precipitates to the ground. In such a case the power D becomes 

D = p cs W ts gH P = P s gH P . (3) 

Here p cs and W ts are respectively the condensate density and the terminal velocity of hy- 
drometeors near the ground surface, precipitation path length Hp is the mean height from 
which the hydro meteors are falling, P s = p cs Wt s is precipitation measured at the surface in 
kg H 2 m -2 s -1 . When a hydrometeor falls to the ground from height Hp, the potential 
energy that it has lost per unit mass is gHp. Precipitation P s tells us how much water 
hits the ground per unit surface area per unit time, so D = P s gHp gives the total rate of 
potential energy loss by all hydrometeors in the column. 

At W t ^> \w\ there is no upward transport of condensate that originated in the lower 
atmospheric layers: the hydrometeors fall to the ground from where they were formed. Rate 
of condensation S (mol H 2 m -3 s _1 ) in the ascending air is equal to 

f dN v dNx ^ _N V p v 



where N„ and N are the molar densities of water vapor and air (IMakarieva et al.l 12010 



Gorshkov et al 



20121 ). The value of S differs from zero only in a certain area z\ ^ z ^ Zz, 
where the relative humidity is close to unity and the water vapor is saturated. Precipitation 
path length Hp is then equal to the mean height of condensation: 

OO Z2 

J S(z)zdz l*(zi)z 1 - 7*^2)^2 + / l*(z)dz 

"^fcr r(*)-r(*l ' (5) 

o 

Here Z\ and z 2 are the heights of the lower and upper boundaries of the condensation area 
(S = at z < z\ and z > z 2 ), 7 = 7* = p* v /p is equal to the ratio of saturated water 



vapor partial pressure p* to air pressure p. (Note that 7* = (M d /M v )q* /[l + q*(M d /M v )] ~ 
(M4/ 'M v )q* , where q* = p* v / pd is the saturated water vapor mixing ratio, M v and M d are the 
molar masses of the water vapor and dry air, respectively. Thus, replacing 7* by q* in ([5]) will 
not significantly affect the estimate of Hp at g* < 1.) Condensation rate 5* normalized by the 
integral in the denominator of ([5]) is the probability that a given hydrometeor reaching the 
surface has fallen from a height between z and z + dz. The last expression in ([51) is obtained 
by integrating the first expression by parts and taking into account that the upward flux of 
air wN is approximately independent of z, that is, d(Nw)/dz = (see Appendix for details). 

We note that Eqs. (j4j) and ([5]) assume that 7 decreases with height solely because of 
condensation that occurs in the rising air parcel; possible mixing with ambient air that can 
lead to a change in 7 is neglected. In the real atmosphere the integration in ([5]) should 
be made over those height intervals where the mean ambient relative humidity is suffi- 
ciently high for condensa tion to occur, e.g., a threshold of Rh(z) > 80% can be applied 



(ILord and Franklin 



1990) • Under the assumption of constant Nw, the vertical profile of 
precipitation P(z) = f 2 S(z')dz' coincides in form with the vertical profile of 7* (z) in the 
area of condensation: we should observe P{z)/P s = r y*(z)/j*(zi) at Z\ ^ z ^ z 2 - Where 
the relative humidity is low but re-evaporation can be neglected, we should observe P(z)/P s 
stays approximately constant independent of changes in j(z). 



3. Numerical estimates of Hp 

The advantage of Eq. ([5]) is that Hp can be estimated from theory. Water vapor condenses 
as the moist saturated air ascends and cools. We allow for the incompleteness of condensation 



which we can define as ( = (z 2 ) / 7* ( z i) > which describes the share of water vapor that has 
not condensed when the saturated air parcel rising from z\ has reached z 2 . For a fully 
saturated atmosphere with surface relative humidity of Rh = 100% we have z\ = 0, z 2 = oo 
in ([5]) and ( = 0. When the vertical distribution of water vapor follows the moist adiabat, 
Hp is unambiguously determined by the surface temperature (see Appendix). In Fig. [IA we 
show the dependence of moist adiabatic Hp on surface temperature T s at a relative humidity 
of R H = 100% and three values of (: 0, 1/2 and 2/3. 

We can see that Hp grows with increasing surface temperature T s , e.g. at ( = we have 
Hp = 3.5 km at 290 K and 5.3 km at 300 K, i.e. Hp increases by around 50% for a 10 
degrees rise in surface temperature. We also find that Hp decreases sharply with increasing 
(: at 300 K and £ = 1/2 we have Hp = 2.4 km, Fig. [TA, i.e. Hp decreases more than 
two-fold compared to the case of complete condensation ( = 0. Height z 2 at ( = 1/2 and 
T s = 300 K is equal to 4.9 km, Fig. [JA. 

In the real atmosphere the lower layer z < Z\ is generally undersaturated with a global 
mean relative humidity at the surface of about Rh = 80%. Value of z\ depends on the 
temperature lapse rate in the lower atmosphere Ti = —dT/dz at z < Z\ (see Appendix). In 
Figs. [IB, HP, and [Tp values of Z\ and a moist adiabatic Hp at ( = are given for R^ of 
80, 60 and 40%, respectively, as dependent on surface temperature for three representative 
values of T^. 5, 6.5 and 9.8 K km -1 . 

We find that the existence of the undersaturated layer z < Z\ does not significantly 
change Hp as compared to the case of Rh = 100%. For example, at T s = 300 K we have 
zi = km and H P = 5.3 km for R H = 100%, Fig. [TA, while for R H = 80% we have 
Z\ = 1.2 km, Hp = 5.6 km at Ti = 5.0 K km -1 and Z\ = 0.46 km, Hp = 5.0 km at 
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r\ = 9.8 K km 1 , Fig. [1)3. The value of z\ grows with decreasing relative humidity. At 
R H = 40%, T s = 300 K and I\ = 5 K km" 1 we have z x = 4.4 km and H P = 7.3 km, Fig. Hp. 
However, under conditions of low relative humidity, the temperature lapse rate cannot be 
much smaller than the dry adiabatic lapse rate, so the estimate of Hp = 7.3 km is not 
realistic. In the realistic case of Fi = 9.8 K km -1 at Rh = 40% we have Z\ = 1.8 km and 
Hp = 4.8 km. In other words, despite the condensation zone being elevated by z\ when 
compared to the case of 100% relative humidity, the mean height of condensation Hp does 
not rise by the same magnitude. This is caused by the temperature dependence of Hp: at 
z\ ~ 1 km, temperature T\ at Z\ is several degrees lower than at the surface, T\ < T s , 
such that the zone of intense condensation is compressed into a smaller vertical space than 
it would be at T\ =T S . So the resulting Hp can be even lower than in the case of z\ = 0. 

It is now possible to evaluate W t from (j2J) at stated values of Hp and P s to see if our 
assumption W t ^> \w\ is realistic. Equating (j2J) and Q we obtain W t = P s Hp/C, where 
C = J °° p c dz is the amount of condensate and W t is the mean vertical velocity of condensate 
in the column. We take the mean precipitation in the tr opical region between 3 0° S a nd 



30° N to be P s = 1.3 m year 1 according to the data of 



Legates and Willmott 



(I1990h . a 



conservative (i.e. low) value of Hp = 3.8 km (this corresponds to T s = 300 K, Rh = 80% 



2002 1. 



and r = 6.5 K km -1 at z > 0, curve 7 in Fig. IB) and C ~ lO^ 1 kg m~ 2 (IWood et al. 
Using these values we obtain W t ~ 1.6 m s -1 . This value is about two orders of magnitude 
larger than typic al time-av eraged large-scale vertical air velocities w < 1 cm s _1 observed in 



the tropics (e.g., 



Rex 



19581 ). We conclude that our assumption Wt \ w\ is reasonable. 



We now consider the case of small terminal velocities, for which Eq. ([3]) does not hold. 



Terminal velocity depends on the size of hydrometeors and turns to zero when the condensate 



particles become vanishingly small. The limiting case is the so-called "reversible adiabat", 
which corresponds to W t = when all condensate travels together with the air and fully 
evaporates in the region where the moist air descends. Note that surface precipitation in 
this case is zero so Eq. ([3]) is inapplicable. 

In order to travel with the air and to reach heights that are significantly larger than 
Hp, the condensate must have a vertical velocity comparable to that of air, w c ~ w and 
Wt <SC \w\. For C ~ 10" 1 kg m -2 and W t <C < 10~ 2 m s _1 in the tropics we obtain 
from (j2j) D < 10~ 2 W m -2 . In other words, even if all condensate in the tropical atmosphere 
consisted of the smallest condensate particles, their contribution to dissipation rate would not 
have exceeded 10~ 2 W m -2 . This means that in the tropics the small amount of condensed 
water that is brought by air updrafts to large altitudes significantly exceeding Hp (JHJ) makes 
little contribution to dissipation rate D ([3]), (j^J), the latter being of the order of 1 W m -2 . 



4. Condensation incompleteness and precipitation effi- 
ciency 

As saturated moist air rises it can mix with drier air from the surroundings. This means 
that its 7 (j3J) drops not because of condensation, but because of dilution through turbulent 
mixing: water vapor is replaced by dry air in the updraft. In such cases condensation is in- 
complete: the water vapor removed from the ascending air by turbulence has not condensed. 
This will affect the value of the so-called precipitation efficiency e. In empirical studies this 
measure, e, is defined as the ratio of precipitation P s at the ground to the inflow of moisture 
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into the updraft. It is commonly observed to be in the vicinity of 20-40% (lFankhauserlll988l ) . 



A common interpretation of these low e value s from o 



Dservation is that it reflects maior re- 



evaporation of condensed water in downdrafts ( iNewton 



1966 



Foote and Fankhauser 



1973|). 



This presumes that all water vapor that has flown into the updraft will condense. Since 
condensation occurs when the relative humidity is equal to unity, this logic would imply 
that the relative humidity within the updraft should remain high up to the layer where 
the water vapor mixing ratio has dropped to a negligible value compared to its value at the 
surface. For example, for a moist adiabat at T s = 300 K a hundredfold reduction corresponds 
to a height of about 14 km, Fig. [2j\. In reality, however, relative humidity drops abruptly 



much earlier - for example, in hurricanes and their a mbient environment it decreases s 



from over 80% to 50-60% at a height of about 4-5 km (ISheets 



1969 



Lord and Franklin 



rarp 



19901) 



y 



The updraft of air carrying the smallest droplets high to the troposphere can continue beyond 
that height, but low relative humidity means that intense condensation cannot. 

In other words, precipitation efficiencies e < 1 do not necessarily imply re-evaporation 
of condensed moisture: rather, a low ratio of precipitation to the water vapor influx can 
indicate an incomplete condensation. Ignoring evaporation, from a simple mass balance 
consideration we have (1 — ()p vs w s = P s , where p vs is the density of water vapor and w s is 
the vertical air velocity at the cloud base, ( = r y*(z 2 )/j*(zi) ~ q* (z 2 ) / q* (^i) is equal to the 
ratio of the water vapor mixing ratio at height z 2 where condensation discontinues and its 
value at the cloud base z = Z\. The flux of water vapor flowing into the cloud is given by 
p vs w s . In this case precipitation efficiency e is equal to e = P s /(p V sW s ) = 1 — C If £ = 1/3, 
this means ( = 2/3, i.e. the condensation zone reaches upward to a height where the water 
vapor mixing ratio decreases by one third as compared to its value at the cloud base. As we 



10 



show in Fig. high values of ( and, hence, low precipitation efficiencies e are associated 
with relatively low precipitation path length Hp. 



5. Comparison with the results of IPauluis et al.l ( l2000l ) 



literature ( IRenno 



2001 



Pauluis et al.l ( 2000 ) based their estimat e of D, the first of this kind in the meteorological 



Pauluis et al. 



20011 ) . on Eq. fl3]). They noted that, if re-evaporation 



is neglected, Hp is equal to the average height where condensation occurs. This is correct, 
b ut some of the s ubseq uent assumptions and derivations are less well justified. Equation (5) 



of 



Pauluis et al. 



( 120001 ) meant to define Hp as Hp = J °° q*dz, where q is water vapor mixing 



ratio, misses the normalization factor q s = q(0) in the denominator (cf. our Eq. (E3) above at 



z i = 0, z 2 = oo and z 2 Y (z 2 ) = 0). Turning to quantitative estimates, 



Pauluis et al. 



( I2000h 



proposed that th e scale height of the saturated water vapor mixing ratio in the tropics is 



about 2.5-3 km. 
refer to 



Pauluis et a 



Emanuel and Bisterl ( 119961 ). Our reading of 



(120001) do not offer a ny clear reasoning fo r this figure but 



Emanuel and Bisterl (1996 



p. 3284) 



finds only a mention of a certain "scale height for water vapor" of around 3 km. However, 
the scale height of saturated water vapor and the scale height of its mixing ratio are different 
atmospheric characteristics. They depend differently on temperature and, hence, height. In 



the tropical troposphere, for example, they range from km to 5 k m and to 10 km, 



'or the 



scale height of water vapor and its mixing ratio, respectively (e.g., 



Fig. 1). 



Makarieva et al. 



2010 



Leading on from their initial suggestions, 



Pauluis et al 



(120001 ) offered several arguments 



as to why Hp in fl3]) should be several times higher than the scale height of the saturated 
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water vapor mixing ratio and takes a value of 5-10 km. First, they noted that the real 
precipitation path length Hp is greater than condensation height obtained from a moist 
adiabat because there is an undersaturated region in the subcloud layer. As we discussed in 
Section 3, while real, the e ffect i s small and does not necessarily lead to an increase in Hp. 



Second, 



Pauluis et al. 



(120001 ) proposed that an increase in Hp can be induced by the 



entrainment of the unsaturated air parcels into the region of saturated ascent. They did 
not, however, specify a mechanism for this. The entrainment of dry air causes the temper- 
ature lapse rate to rise above the moist adiabatic value, such that the temperature drops 
more rapidly with height. If, despite the dry air entrainment, the condensation nevertheless 
continues in the updraft prompted by this additional cooling in the ambient environment, 
the change in local lapse rate will reduce rather than raise the mean condensation height 
Hp. Curve 7 in Fig. IB illustrates the dependence of Hp on T s for Rh = 80% and a mean 
tropospheric lapse rate of 6.5 K km -1 instead of moist adiabatic lapse rate. In this case Hp 
depends little on surface temperature and ranges between 3 and 4 km. If, on the other hand, 
condensation is discontinued by the drop in relative humidity associated with the removal 
of water vapor and its replacement by dry air in the region of ascent, then Hp is limited by 



the height where t 



Pauluis et al 



re dr y air entrainment occurred. In neither case does Hp increase. 



( 120001 ) also mention that the real precipitation path length is increased 



by the fact that some hydrometeors are lifted by updrafts to high altitudes. However, as 



discussed Section [21 such hydrometeors are small and po ssess terminal 



much smaller than air velocity, such that Eq. ((2]), which 



Pauluis et al. 



velocities that are 



(120001 ) intended to 



use, is thus inapplicable. On average, these hydrometeors make only a negligible contribution 
to the total dissipation rate associated with precipitation owing to both their slow terminal 
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velocity and their small combined mass. 



The final argument put forward by 



Pauluis et al. 



(J2000|), and the only quantitative one, 



concerns re-evaporation. They presume that t his e ffect can lead to a significant u n deres- 



tim ate of the rea l value of Hp. 
and 



Pauluis et al. 



Ferrier et al. 



( 120001 ) quote the work of 



Fankhauserl dl988h 



(119961 ) to support the statement that a significant part (from half to two- 



thir ds) of all condensed m oisture actually re-evaporates and does not hit the ground. From 



this 



Pauluis et al 



( 120001 ) suggest that if evaporation occurs uniformly as the hydromete- 



ors are falling, this process increas es the effect i ve pre cipitation path length by a factor of 



1.5-2. We note for the record that 



Fankhauserl (119881 ). who investigated empirical data on 



the water budget of convective clouds including precipitation efficienc y, does not mention 



any quantitative estimate of the re-evaporation of condensed moisture. iFerrier et al.l (11996 . 
p. 2105), on the other hand, do report the magnitude of re-evaporation as compared to total 
condensation within a squall, but their results come from a numerical model rather than 
observational evidence. In order to estimate the actual rate of evaporation of condensate in 
the downdrafts we would need to perform careful estimates of condensate transport within 
the cloud - rather than measuring the transport of total moisture that is dominated by water 
vapor. Estimates of sufficient accuracy are not available. We can, in contrast, be confident 
that the effect of incomplete condensation associated with low precipitation efficiency is real. 
But, as discussed in the previous section, this will decrease rather than increase the estimate 



oiH 



Furthermore, even if evaporation in downdraft s did constitute a significant fraction of 



total condensation, the suggestion of 



Pauluis et al. 



( 120001 ) about a 1.5-2 increase in effective 



Hp would still be incorrect. Let us first see how this conclusion was reached, because 
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Pauluis et all ( 120001 ) do not explain this in any detail. The argument of iPauluis et al.l ( 120001 ) 
derives from Eq. fl3]) along with the one-dimensional continuity equation for condensate 
particles. If jp = p c {z)Wt is the downward flux of condensate at po int z, then the contin uity 



Pauluis et al. 



(120001 ). we 



equation is djp jdz — E, where E > is evaporation. If, following 
assume that E is constant and that jp(0) = (1/3) jp(Hp) (evaporation has decreased the 
original precipitation flux by two thirds as it traveled from z = Hp to z — 0), then we have 
jp(z) = jp(0)(l + 2z/Hp). This allows us to calculate D from (J2J) as D = J^ p jpqd 



2jp(0) qHp = 2P s gHp. Had this logic been correct, we could conclude, as did 



Pauluis et al 



(120001 ) . that Eq. ([3]) indeed underestimates the real dissipation by half. 



There are, however, two errors in this reasoning. The first one consists in neglecting 
the fact that local dissipation rate p c gWt and local evaporation with respect to droplet 
size behave differently. Since absolute evaporation rate is proportional to droplet area, 
the smallest droplets evaporate most rapidly. If we have equal amounts M (g) of small 
droplets with radius t\ and large droplets with radius r 2 > r±, the rate of depletion of 
total condensate from the evaporation of small droplets will be rijr\ times faster than from 
large droplets (evaporation rate oc Ns oc (M/m)r 2 oc M/r, where s oc r 2 is droplet's 
surface area, N oc M/m oc M/r 3 is the number of droplets, m oc r 3 is droplet mass). In 
comparison, because of the fact that terminal velocity grows with increasing droplet radius, 
the contribution of the small droplets to total dissipation will be lower than that of the large 
droplets. In theory, for spherical droplets with Wt oc r 2 , it will be lower by (r 2 /ri) 2 times. 
For example, with r 2 /ri = 10, 90% of all evaporation will come from droplets that make a 
1% contribution to total dissipation. 

The second error in their reasoning is that they implicitly use jp = p c (z)Wt to rep- 
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resent the downward flux of condensate. This neglects the important role of vertical air 



movements in transporting the smallest condensate particles. Rather they should have used 
jc = Pc{z){w — Wt), where w < is the downward velocity of air. The w term is particularly 
important in consideration of evaporation, because the smallest droplets with W t <C —w are 
so slow that they can be only transported by the do wndraft. This means that the continuity 



Pauluis et al. 



fbOQOl ) is not valid. The 



equation djp/dz = E underlying the reasoning of 
correct equation dj c /dz = E does not offer any insights regarding the frictional dissipation 
of D in the column because the distribution of w remains unknown. 



To summarize, we find no support for the claim of 



Pauluis et al. 



(120001 ) that precipitation 



path length should be several times higher than the value of Hp given by 



6. Numerical estimate of D 

We have shown that Hp grows with increasing incompleteness of condensation (, increas- 
ing surface temperature and decreasing lapse rate I\ at z < z±, Fig. [TJ Among these, the 
incompleteness of condensation ( is both the least known and the most influential, Fig. 1A. 
It is closely linked to convection depth. We suggest that the uncertainty of the mean D 
values for the tropical region is largely determined by the uncertainty in Hp, which is, in 
its turn, largely reflects uncertainty of (. Taking T s = 300 K as the mean surface tempera- 
ture during precipitation in the tropical region, Ti = 5 K km -1 and Rh = 80%, we obtain 
Hp = 5.6 km for complete condensation ( = 0, Fig. dJ3. At the lower end is the estimate of 
Hp obtained assuming ( = 2/3 (corresponding to a low precipitation efficiency of e — 1/3), 
which at T s = 300 K, T 1 = 5 K km" 1 and R H = 80% gives H P = 2.5 km. 

15 



Pauluis et al.l (120001 ) estimated D from the mean latent heat flux Q instead of precipita- 
tion P = Q/L v in the tropics considering that D/Q = gHp/L v , where L v = 2.5 x 10 6 J kg -1 
is the heat of vaporization. Using Q = 100 W m -2 and Hp values of 2.5 and 5.6 km we 
obtain a range of 1.0-2.2 W m -2 for the mean tropical value of D. 

For the global mean temperature T s = 288 K at Rh = 80% and T\ = 6.5 K km^ 1 we 
have Hp = 3.6 km for £ = 0, Fig. [TJ3, and Hp = 1.5 km for C = 2/3. U sing the mean value 



of 2.5 km and mean global precipitation of 1 m year 1 (ILVovitch 



19791 ) we obtain a global 



mean value of D ~ 0.78 W m~ 2 from This means 4 x 10 14 W for Earth as a whole 
and 1.2 x 10 14 W for the gravitational power of precipitation on land. (We note that the 
estimate of 10 14 W for land was prev iously obtained based on Eq. ([3]) in a different context 



discussing renewable energy sources (IGorshkov 



1982 



p. 6). 



If the vertical profile of precipitation P(z) is known, the value of D can be estimated 
directly from (121) under the assumption that P(z) = p c (z)W t . This was recently done by 



Pauluis and Diasl (120121 ) who used satellite-derived P(z) profiles from the Tropical Rainfall 



Measurement Mission and es timated D for the t ropica l region between 30° S and 30° N to be 



Pauluis and Diasl (120121 ) went on to estimate Hp from Eq. 



1.8 W m 2 . Using this value 
by dividing D (J2j) by P s q, Hp = D/(P s q). H aving obtained values of 5.1 km for the ocean 



and 6.9 km for land, 



Pauluis and Diasl (120121 ) con cluded that these re sults agree with their 



earlier estimates of 5-10 km for Hp in the tropics (IPauluis et al. 



2000). They proposed that 



the larger value obtained over land indicates more intense convection than over the ocean. 

The derivation of precipitation rates from the satellite radar data involves considerable 
uncertainties. The estimation process involves a large number of empirical relationships be- 
tween reflectivity and precipitation rates as well as various assumptions concerning the prop- 
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erties of the hydrometeors like their size distribution and termin al velocity (e.g., lUijlenhoet 



2001 



Durden et al. 



1998 



Bowman 



2005 



Prat and Barros 



rithms perform differently over land than they do over oceans (IPrat and Barros 



20091). The common 



also perform differently on the ground surface versus at the top of clouds (IDurden et al. 



y use d algo- 



200 91). They 



1998|) • The estimate of 1.8 W m -2 reported without any quantified assessment of the as- 
sociated uncertainties ranges does not constrain D any more accurately than do theoretical 
estimates, although it ca n be noted that this e stimate falls out of the 2-4 W m -2 established 



for the tropical region by 



Pauluis et al 



fl200Q[ ). 



Moreover, our examination of the precipitation profiles P{z) shown in Fig. 2 of 



Pauluis and Dias 



(120121 ) raises further problems. These data support neither the estimate of D = 1.8 W m~ 2 
reported for the tropical region as a whole, nor the estimates of Hp made for land and ocean. 
Indeed, integrating these profiles yields a value of 1.4 W m -2 for all the three profiles (the 
tropics as a whole, land and ocean). This 30% discrepanc y raises further q u estion s con- 



Pauluis and Dias! (120121 ). The 



cerning the validity of the numerical estimates reported by 
discrepancy illustrates the potential inaccuracy associated with using satellite radar data on 
precipitation. Heights Hp estimated from the corrected values of D using P s = 2.6 mm day -1 
for ocean and P s = 2.1 mm day -1 for land yields Hp = 4.7 km for the ocean and Hp = 5.7 km 



for land instead of, respectively, 5.1 km and 6.9 km obtained by 



Pauluis and Diasl fj2012h . 



Furthermore, the derivation of Hp from Eq. Q is potentially misleading, as such an 
estimate is sensitive to the estimate of surface precipitation. Satellite radar based assessments 
of precipitation are based on the reflectivity coefficient of precipitating particles. The radar 
is unable to accurately distinguish surface precipitation because of the high reflectivity of 
the surface. Note for example, that if surface precipitation is underestimated in the lowest 
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kilometer, this may make little impact on the column-integrated D (j2J), but will have a 
considerable impact on the value of Hp. (We also note that precipitation rate s in the upper 
part of the atmosphere can be, on the contrary, overest imated by the radar (IDurden et al. 



1998 ).) According to Fig. 2 of 



Pauluis and Diasl (120121 ). precipitation in the lowest 1 km 



makes about a 1/5 contribution to the column-inte grated value o. 



D. If surface precipit ation 



1998 



Bowman 



20051 ) and 



in this lowest region is underestimated by 25% (IDurden et al. 
constitutes 75% of the real value, this corresponds to a 25%/5 = 5% underestimate in total D. 
But Hp = D/(gP s ) is then overestimated by a factor of (100 — 5)%/75% = 1.3. Applying this 
additional correction factor to our corrected Hp estimates we obtain Hp = 4.7/1.3 = 3.6 km 



instead of 5.1 km for the ocean and Hp = 5 .7/1.3 = 4.4 



on instead of 6.9 km for land. 



Both values are outside the 5-10 km range of 



Pauluis et al. 



(12000l ). but within the 3-5 km 



range estimated from our theoretical analysis. We also note that these values are close to 
Hp = 3.8 km that at 80% surface relative humidity is characteristic of the mean tropospheric 
lapse rate 6.5 K km -1 in a saturated atmosphere at T s = 300 K, see curve 7 in Fig. [IB- 



Tropical land include some exceptionally wet regions like the Am azon and Congo 



where precipitation is 2-3 times higher than over the nearby ocean (IMakarieva et al. 



orests, 



2012 ). 



But the tropics also include very dry regions such as the Sahara desert and the Australian in- 
terior. Combining these wet and dry regions under one and the same catego ry and calculating 



a sing le mean precipitation profile for all tropical land, as done in Fig. 2 of 



Pauluis and Dias 



(120121 ). is of questionable value given the diversity of physical settings. (Note that these 



concerns have less significance for a tropic-wide estimate the values of which are dominated 
by the oceans). In the driest regions of the Earth where surface precipitation tends to 



zero, estimating Hp from surface precipitation lacks any physical meaning, as Hp — > oo 



at P s — > 0. Thus any estimated value for Hp does not carry any information about the 



real vertical distri 



jution or intensity of precipitation. We note in this context that while 



Pauluis and Diasl ( 120121 ) concluded that high values of Hp over land are indicative of a more 



intense convection, judging from their Fig. 3 one might also think that the region of most 
intense convection on Earth is the inner part of Sahara desert. Here Hp, as estimated with 



use o f the near zero value of surface precipitation, reaches beyond 10 km ( iPauluis and Dias 



2012 



Fig. 3). This apparently nonsensical result illustrates the need to have a consistent 
theoretical basis for any analysis of empirical evidence. 

As we discussed in Section |2j in a saturated atmosphere the vertical profile of precipitation 



P(z)/P s should coincide with the vertical profile of j*(z)/j(0). I n Figs. [2)3, [2]C, an d [2p 



Pauluis and Diasl (120121 ) is 



the 



mean vertical profile of tropical precipitation taken from Fig. 2 of 
contrasted against theoretical profiles of 7*(^)/7(0) calculated for different values of surface 
temperature, relative humidity and temperature lapse rate. We see that saturated moist 
adiabatic profiles of 7*(^)/7(0) tend to overestimate P(z)/P s in the upper atmosphere, 
Figs. 03 and|2p. The observed mean profile P(z)/P s is confined between 7*(z)/7(0) profiles 
built for surface relative humidity from 80% to 40% and having a constant lapse rate of 
T = 6.5 K km" 1 , Fig. 03. 



7. Discussion 

In order to investigate how an atmospheric phenomenon responds to changes in atmo- 
spheric parameters it is important to establish a sound theoretical basis concerning the key 
physical relationships. In this paper, building from basic physical principles and relation- 

19 



ships, we evaluated the rate of turbulent frictional dissipation associated with precipitation. 
We discussed how precipitation path length, Hp, is the key parameter controlling this rate, 
and investigated how it depends on surface temperature, humidity and the vertical extent 
of the area where condensation occurs. 

We now consider how the frictional dissipation relates to the dynamic power of atmo- 
spheric circulation. We can illustrate this relationship with a simple example. Consider a 
hanging weight tied with a short rope to an extended spring. The system is in a state of 
equilibrium with the weight exactly balanced by the tension of the spring. When we cut 
the rope the weight can fall. At this moment two types of potential energy have become 
available: the first is the weight's energy in the gravitational field and the second is the 
energy of the stretched spring. The first potential energy is transformed to kinetic energy 
as the weight falls, while the second potential energy is transformed to kinetic energy as the 
spring accelerates upward. It is clear from this example that the two potential energies are 
independent in magnitude: the first depends on the initial height of the weight above the 
surface, the second depends on the elasticity of the spring and is independent of height and 
of what happens to the weight after the rope is cut. 

Likewise upon condensation, two types of potential energy are formed: first the potential 
energy of falling hydrometeors as dictated by the precipitation path length Hp (JSJ) and 
second the potential energy of the non-equilibrium pressure gradient that results from the 
disappearance of water vapor from the gas phase. The key peculiarity of this second process 
consists in the fact that this potential energy is coupled to the vertical motion of moist air 
and is released only when the air moves upwards. 

Since precipitation and condensation always accompany each other, it is of interest to 
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compare the powers associated with each of the two processes and how these might influence 
atmospheric motion. Previously we have argued that the dynamic power of condensation- 
induced circulation D c (W m -2 ) per unit surface area can be estimated as D c = P S RT/M V , 
where R = 8.3 J mol" 1 K -1 is the universal gas c onstant and T is the mean temperature in the 



atmo s pheric column where c ondensation occurs (IMakarieva et al 



2011 



Gorshkov et al. 



2010 



Makarieva and Gorshkov 



20121 ). This formula results from the proposition that the dynamic 
power T> c (W m -3 ) per unit air volume of the upward pressure gradient force induced by 
condensation and associated with the non-equilibrium vertical gradient of saturated water 



vapor, is equal to 



w 



dp* v dp 
dz 



^ dz 



wp- 



<9 7 * 
dz 



RTS, 



(6) 



where S = wNd'-f* /dz is the net condensation rate per unit volume in mol H2O m~ 3 s" 1 , N 
(mol m -3 ) is molar density of air, see 01]). Neglecting the minor dependence of T ~ T s on z 
and observing that the integral of S over z is equal to P s /M v , see Appendix, we have 

P,RT, 



D, 



T> r dz 



M„ 



To illustrate the value of this theoretical result we can now estimate the globa 



cond e nsation-induced circulation. Using a global mean value of P s of 1 m year 1 (jL'vovitch 



1979 



Legates and Willmott 



(7) 



power of 



19901 ) and global mean surface temperature of T s = 288 K, 
we have D c = 4 W m -2 . For the tropical region with mean P s = 1.3 m year -1 we have 
D c = 5.5 W m~ 2 . Since under conditions of hydrostatic equilibrium the work done by the 
vertical pressure gradient is compensated by the work done by gravity, the kinetic energy of 
the large-scale air flow derives from the horizontal pressure gradient alone. The power of this 
horizontal force per unit air volume is equal to udp/dx, where u is the horizontal velocity 
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component parallel to the horizontal pressure gradient. The dynamic power of atmospheric 
circulation (the rate at which the kinetic energy is generated) can therefore be estimated 
from the observed horizontal pressure gradients and the observed u values. It can also be 
estimated as the power of turbulent dissipation of the air flow under the assumption that 
in the stationary case the dynamic power that creates the kinetic energy is equal to the 



power of turbulent dissipation D t of this energy. The available globa 



based estimates o 



Peixoto and Port 



thes e powers are in the range of 2-4 W m 2 (jOort 



mean observation- 



1964 



Lorenz 



1967 



1992). Our global mean estimate, derived from basic principles, falls at 
the upper edge of this range. 

The ratio of the powers of the two processes - the frictional dissipation power D of 
hydrometeors and the dynamic power D c of condensation-induced circulation - is given by 

D _H P fc RT S 

D c ~ h v > hv ~ M v g W 

This ratio does not depend on precipitation rate P s but grows with temperature owing to 
the dependence of Hp on T s , Fig. [TJ The value of h v OH]) has the meaning of the scale height 
of (unsaturated) water vapor in hydrostatic equilibrium, at T s = 300 K it is equal to 14 km. 
With the large-scale values of Hp not exceeding 6 km at mean surface temperatures not 
exceeding 300 K, Fig. [TJ we obtain a general estimate of D/D c < 0.4. 

Noting that D grows with surface temperature, Fig. 1, we can estimate when D exceeds 
D c . We do this using (jSJ) and assuming that the budget of energy turnover for a condensation- 
induced circulation has the form of D c = D t + D, where D t is turbulent dissipation of the 
large-scale air flow. At R H = 80% and Ti = 5 K km^ 1 D equals D c at T s = 323 K, Fig. [IB. 
At higher temperatures in a moist adiabatic atmosphere any significant circulation due to 
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condensation will be prevented because of the insufficient dynamic power to overcome the 
energy losses associated with frictional dissipation due to precipitation. If the atmosphere is 
not moist adiabatic but has a constant lapse rate of T = 6.5 K km -1 , D grows much more 
slowly with increasing surface temperature (see Fig. 1, curves 7). It does not approach D c 
anywhere at T s < 360 K, i.e. in the entire range where the approximation 7 < 1 on which 
Eq. ([5]) is based holds. From these considerations we conclude that frictional dissipation 
due to precipitation is insufficient to arrest condensation-induced atmospheric circulation on 
Earth. 
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APPENDIX 



Equations for calculating z\ and Hp 



Eq. d5J) is obtained using (jlj) and considering that 



z-1 



S(z)zdz = wN^*z 



wN-f*dz+ / — '-rfzdz 



z-i 



wN < 7*2; 



z\ 



Z-2 



+ / J*dz 



S{z)dz = wNY 



zi 



Z2 



d(wN) 
dz 



Zl 



■y*dz 



dz 



(Al) 



(A2) 



From the one-dimensional stationary continuity equation we have d(wN)/dz = —S. This 
means that the terms discarded in (1A1I) and (1A2I) constitute a small magnitude of the order 
of 7* 1 as compared to the initial terms. Quantity wN changes little with z (by a relative 
magnitude of the order of 7*) as compared to 7* that changes severalfold. Therefore Nw 
can be assumed to be constant and cancelled from both the denominator and nominator in 
ratio (I31). The inaccuracy of the resulting expression for Hp is of the order of 7* and, for 
temperatures of interest, does not exceed 10%. 

The system of equations for moist adiabat solved to calculate Hp in Fig. [TA is as follows 
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flMakarieva et all 120 Id iGorshkov et all 12012 ): 




(A3) 



(A4) 



1 dp Mg 



(A5) 



pdz RT 



where temperature T, air pressure p, and 7* =p*/p- the relative partial pressure of saturated 
water vapor - are functions of height z. Here £ = L v j (RT), L v = 45 x 10 3 J mol" 1 is heat of 
vaporization, R = 8.3 J mol -1 K" 1 is the universal gas constant, M = (1 — 7*)M^ + 7*M„, 
Ma = 29 g mol ~ 1 , M v = 18 g mol" 1 , = i?/c p = 2/7, c p is the molar heat capacity of air at 
constant pressure (J mol -1 K" 1 ). Equation ( 1A3|) results from the first law of thermodynamics 
for moist air saturated with water vapor. Equation flA4[) derives from the definition of 7* 
combined with the Clausius-Clapeyron law. Equation (1A5I) is equivalent to the condition of 
hydrostatic equilibrium for ideal gas. 

The boundary conditions for the surface z = at a given surface temperature T s read 



where p* Q and £0 = L V /(RT ) correspond to some reference temperature T . We take T = 
303 K, 

PvO ~ hPa and the standard value for the atmospheric pressure p s — 1013 hPa. 
The dependence of vaporization heat L v on temperature is neglected. In this case £0 = 18. 

For a fully saturated atmosphere Z\ — in ([5]). Numerical evaluation of the system of 
Eqs. (IA3j) - flA5[) allows us to obtain the unknown functions T(z), p(z), 7* (2) and to calculate 



T = T. 



(A6) 



P = Ps, 



(A7) 



P*v( T ) =P*vO eX V(£o -0. 



(A8) 
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Hp (jSJ) as shown in Fig. [QA.. 

In Figs. HB, CEP, and [Tp the atmosphere at the surface z = is not saturated and has a 
relative humidity of 80%, 60% and 40%, respectively. To find height z\ where the relative 
humidity reaches unity, we assume that within the range ^ z ^ Z\ the non-saturated 
7 = Pv/p is constant and temperature T(z) drops versus height 2 at a constant lapse rate 
Ti, T(z) = T s — Y\z. Then the non-saturated pressure p v of water vapor is given by 

p v (z)=pl(T s )R H ^, (A9) 

Ps 

where Rh is relative humidity at the surface z — 0, saturated pressure p* of water vapor is 
governed by the Clausius-Clapeyron law (IA8I) and pressure p(z) conforms to the condition 
of hydrostatic equilibrium ( )A5j) with M ~ M^. Height Z\ where relative humidity becomes 
unity and condensation commences, is a function of the surface temperature T s . We find 
height Z\ as the solution of equation 

p v ( Zl ) = p* v (T x ), T ± = T 8 - T lZl . (A10) 

The atmosphere is assumed to be saturated and moist adiabatic within the range Z\ ^ 
z ^ z 2 . We need to evaluate the system of Eqs. ( 1A3I) - (1A5I) in order to find functions T(z), 
p(z), and 7* (z). However, the boundary conditions should now be imposed not at the ground 
surface but at z — z±, so that T = T\ and p = p(z\) instead of T s (IA6I) and p s (IAT[) at z — 0. 
Using the obtained solutions, condensation rate is calculated from 

Curves 7 in Figs. [IB; HP; and[Tp are obtained for an atmosphere that is unsaturated at 
z < Z\, saturated at z ^ z\ and has a constant temperature lapse rate of V — 6.5 K km -1 
at z ^ 0. The solution at z ^ Z\ is obtained by solving the system of Eqs. (IA4I) . (|A5[) and 
BT/dz = -r instead of flA3|) . 
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List of Figures 

1 Precipitation path length Hp as a function of surface temperature T s at var- 
ious values of surface relative humidity Rh- 100% (A), 80% (B), 60% (C), 
40% (D). 

A: Moist adiabat in a fully saturated atmosphere, z\ = 0. Height of conden- 
sation area z 2 is equal to the height where ( = 7*(z 2 )/7*(^i) = (triangles), 
C = 1/2 (circles) or £ = 2/3 (squares). 

B, C, D: Curves 1-3 denote z±, curves 4-7 denote Hp at £ = 0. Tempera- 
ture lapse rate ri at z < Z\ is equal to 9.8 K km^ 1 (red squares), 6.5 K km^ 1 
(blue circles, stars), 5 K km -1 (green triangles). At z ^ z\ the atmosphere 
is moist adiabatic, except for curves 7, which are for an atmosphere that is 
saturated above %\ and has a constant temperature lapse rate T = 6.5 K km -1 
everywhere at z ^ (see Appendix for details). 
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Vertical distribution of 7*(z)/7(0) under different atmospheric conditions. 
Dashed curve denotes the mea n obse rved tropical precipitation P(z)/P(0) 



from Fig. 2 of 



Pauluis and Diasl (120 12 ). 



A: Moist adiabat in a completely saturated atmosphere, Rh = 100%, cf. 
Fig. 1A. 

B: The atmosphere is unsaturated at z < z\ (Rh = 80%, Ti = 5 K km -1 ) 
and moist adiabatic at z ^ z±, cf. Fig. [IB- 

C: The atmosphere is unsaturated at z < Z\ and moist adiabatic at z ^ z±. 
D: The atmosphere is unsaturated at z < z±, saturated at z ^ Z\ and has a 
constant temperature lapse rate T = 6.5 K km -1 at z ^ (cf. curves 7 in 
Figs. [IB, HC, and HP). 
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Fig. 1. Precipitation path length Hp as a function of surface temperature T s at various 

values of surface relative humidity R H : 100% (A), 80% (B), 60% (C), 40% (D). 

A: Moist adiabat in a fully saturated atmosphere, z\ = 0. Height of condensation area Z2 

is equal to the height where ( = 7*(z 2 )/7*(^i) = (triangles), ( = 1/2 (circles) or £ = 2/3 

(squares). 

B, C, D: Curves 1-3 denote z±, curves 4-7 denote H P at £ = 0. Temperature lapse rate 
r\ at z < Zi is equal to 9.8 K km -1 (red squares), 6.5 K km -1 (blue circles, stars), 5 K km -1 
(green triangles). At z ^ z\ the atmosphere is moist adiabatic, except for curves 7, which 
are for an atmosphere that is saturated above z\ and has a constant temperature lapse rate 
T = 6.5 K km -1 everywhere at z ^ (see Appendix for details). 
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Fig. 2. Vertical distribution of 7*(z)/t(0) under different atmospheric conditions. 
Dashed curve denotes t he mean observed tropical precipitation P(z)/P(0) from Fig. 2 of 
Pauluis and Diasl ( l2012t ). 

A: Moist adiabat in a completely saturated atmosphere, Rh = 100%, cf. Fig. 1A. 

B: The atmosphere is unsaturated at z < z\ [Rh = 80%, I\ = 5 K km -1 ) and moist adia- 

batic at z ^ z±, cf. Fig. [TJ3. 

C: The atmosphere is unsaturated at z < z\ and moist adiabatic at z ^ Z\. 

D: The atmosphere is unsaturated at z < z\, saturated at z ^ Z\ and has a constant 

temperature lapse rate V = 6.5 K km -1 at z ^ (cf. curves 7 in Figs. [1)3, [Tp, and[Tp). 
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